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Selected, recent results, primarily from collider experiments but including some fixed 
target experiments, are presented as a survey of Quantum Chromodynamics (QCD). The 
concepts of leading order and next-to-leading order QCD are introduced. Inclusive pp jet 
and dijet production and deep inelastic ep scattering at very large momentum transfer 
are shown to be in good agreement with pcrturbativc QCD (pQCD). Dijet, three-jet 
and multi-jet results from pp, ep, and ee colliders at moderate Q 2 are also compared to 
pQCD. BFKL searches from all three colliders are discussed. Recent measurements of 
structure functions and contributions to the parton distribution functions are presented. 
New measurements of a s are summarized, the world average is ct s = 0.1184 ± 0.0031. 

1. Introduction and a Brief Introduction to QCD 

Quantum chromodynamics (QCD) beautifully describes the strong interaction over 
a surprisingly broad range of phenomena - from the non-perturbative regime of 
hadronic structure to the highest energies observed in particle colliders. Recent mea- 
surements of high momentum transfer jet production and deep inelastic scattering, 
multi-jet production, proton parton distribution functions, and the strong coupling 
constant convincingly illustrate the breath of experimental progress. These results 
have, in turn, stimulated impressive theoretical progress and sensitive searches for 
new interactions and dynamics. Because of the immense quantity of QCD research 
underway the results discussed here are necessarily and unfortunately incomplete, 
this document might best be viewed as an introductory overview. 

The proton-antiproton interaction, a general scattering process, nicely intro- 
duces the concepts of leading order and next -to-leading order perturbative quan- 
tum chromodynamics. Inelastic scattering between a proton and an antiproton can 
be described as an elastic collision between a single proton constituent and single 
antiproton constituent. These constituents are collectively referred to as partons 
and in QCD are quarks and gluons. Predictions for jet production are given by fold- 
ing experimentally determined parton distribution functions / with quark and gluon 
two-body scattering cross sections a. The two ingredients can be formally combined 
to calculate any cross section of interest: a = I dx\dx2fi{xi, fJ<p)fj(%2, H 2 F ) 
a[x\P, X2P, a s (//|j.), Q 2 1 ' [i 2 F , Q /fijj}- The nonperturbative parton distribution func- 
tion fk(xi,fi 2 7 ) describes the momentum fraction x of the beam momentum P car- 
ried by the I th parton of type k. The hard two-body interaction between the gluon 
and quark partons can be calculated with perturbative QCD (pQCD) and is a func- 
tion of the perturbative or strong coupling constant a s , the hard scale momentum 
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transfer between incoming particles Q, and the renormalization scale 

The factorized scattering has been illustrated with a leading order (LO) process; 
that is, a process with the minimum number of vertices (or coupling constants) to 
describe two final states. Although useful, the leading order picture (where one 
parton results in one jet) is too simple and has an unphysical dependence on 
Next -to-leading order (NLO), or for this two jet process O(af) calculations, include 
additional parton emission and have reduced sensitivity to [Ir. Depending on the 
proximity of the outgoing partons, a "jet" could result from one or the combination 
of two partons. NLO calculations crudely model fragmentation, thereby obviating 
the need for fragmentation functions. 



2. QCD at the Highest Momentum Transfers 



A complete theoretical description of inclusive jet production, pp — *■ j -\-X, requires 
proper treatment of the final state radiation and accurate measurements of the 
parton distribution functions, pdf's. The inclusive jet cross section is reported 
as d 2 7]/ 'dErdrj. Et — EsinO where E is jet energy and 9 the angle between the 
proton direction and the jet. The pseudorapidity, is defined as —ln(tan(0/2)). 
Kincmatically, an individual jet is characterized by Et , ?7, and azimuth 0. Jets are 
found by clustering energy in a cone of radius 0.7 in 77 — <f> space EJ. 

The inclusive jet cross section has been intensively studied at beam energies of 
900 GeV/c at the Fermilab Tevatron proton-antiproton collider in Batavia, Illinois 
by the D0 collaboration for \rj\ < 3.0 i'B and by the CDF collaboration in the 
region 0.1 < \rj\ < 0.7 B@. The data now span nearly ten orders of magnitude for 
jet energies between 20 and 450 GeV. The percentage difference as a function of Et 
between the data and theory normalized to the theory for the D0 data is shown in 
Fig. |l|. Events at the highest Et correspond to x values of 0.5 and Q 2 greater than 
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Figure 1: The difference between the inciusive cross section and NLO theory nor- 
malized to NLO theory. 
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10 5 GeV 2 . Data points include statistical errors only and the bands indicate the 
magnitude of the systematic errors. The figure includes a NLO prediction using 
the CTEQ4HJ pdf. There is excellent agreement at all Et and r\ and, in particular, 
no indication of excess production at large Et ■ The CDF collaboration finds similar 
agreement in the region O.f < \r\\ < 0.7, but reports a 25% discrepancy □ for Et 
above 400 GeV. 

The good agreement is in sharp contrast to an earlier, lower statistics results 
from the CDF collaboration which suggested excess production at large Et relative 
to contemporaneous NLO predictions □. The early CDF results prompted a closer 
look at theoretical and pdf uncertainties contributing to the predictions. In fact, 
inspired by the initial CDF result, the CTEQ4HJ pdf used in Fig. |l| includes a 
strengthened high-x gluon component. All D0 and CDF results are consistent 
within statistical and systematic errors 

The dijet mass distribution for the leading two jets of a pp — > j + j + X event at 
beam energies of 900 GeV also constitutes a sensitive search for new physics. The 
cross section is reported as d 3 i]/ 1 dMjjdr\\dr\i where Mjj is the invariant mass of the 
leading two jets. Both D0 and CDFhaye published results which show good agree- 
ment between data and NLO QCD UBu. As shown in Fig. |[ the DO collaboration 
has taken the ratio of the central dijet mass distribution \rj\ < 0.5 to the forward 
distribution, 0.5 < \rj\ < 1.0. This distribution is sensitive to quark compositeness 
because the associated jet production will be predominantly central in rapidity. In 
a manner completely analogous to Rutherford scattering, excess jet production at 
very large transverse energies signals the presence of quark compositeness. The 
curves in Fig. || include NLO predictions with and without a additional jet produc- 
tion due to compositeness. The solid curve postulates no substructure, and DO has 
set a limit of 2.4 TeV or ~ 10~ 19 m on any substructure 

HERA, an electron-proton collider in Hamburg, Germany provides a comple- 
mentary and comprehensive opportunity to examine QCD over a very large range 
of momentum transfer (0 to 3 x 10 4 GeV 2 ). The introductory description of pp jet 
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Figure 2: Comparisons of central to forward dijet cross section ratios as a function 
of mass to theoretical predictions. See text for details. 
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production can be modified to ep scattering if one incoming hadron is replaced with 
an electron and the exchanged gluon by an electroweak boson . A neutral current 
(NC) reaction is characterized by an exchanged photon or Z, while charged current 
(CC) reactions involve an exchanged W. The CC and NC cross sections are re- 
ported as d 2 a/dQ 2 . HERA ran from 1994-1997 with positrons at a center-of-mass 
energy of 300 GeV and from 1997-1998 with electrons at an energy of 320 GeV. 

As shown in Fig. |L the ep — > e + X neutral current cross section spans eight 
orders of magnitude 119. Below Q 2 — 10 4 GeV 2 the data are well described by 
NLO theory. However, the NC e + p data shows excess production at very high Q 2 , 
an effect observed (and which generated great ^xcitement) by both the HI and 
ZEUS collaborations in the 1994-1997 data set BH . In contrast, the 1997-1998 
erp NC cross_section (also shown in Fig. |^) and similar results from the ZEUS 
collaboration 113 are both well described within statistical errors at all Q 2 by NLO 
QCD. Although not shown, CC ep — > v 4 | _JC_cross sections for both electrons and 
positrons are also well described by QCD I13'EJ. 

The pp and ep colliders both test the limits of QCD and the Standard Model 
(SM) at the highest momentum transfer. In the mid-nineties hints of excess jet and 
NC events at large Et or, equivalently, Q 2 were observed at both colliders. However, 
since then new results (characterized by higher statistics or new projectiles) and 
modified theoretical predictions (characterized by new pdfs), strongly suggest the 
SM can currently describe the data. The lure of new physics is exciting and searches 
for high Q 2 discrepancies will continue, especially with the ever increasing statistics 
at HERA and Fermilab. However, as seen in the next section, the lower Q 2 regime 
of multi-jet production proves an interesting window onto the behaviour of QCD. 

3. Multijet Production 

All three colliders, HERA, the Tevatron, and the LEP e + e~ machine in Geneva, 
Switzerland prove prolific laboratories for the study of multi-jet production. At 
HERA and LEP jets are conventionally found with recombination algorithms which 
cluster energy according to their proximity in Et , r\ , and azimuth El. At LO, HERA 
dijet production occurs through the QCD Compton graph (where a photon from 
the incoming lepton scatters off a quark from the hadron to produce a gluon jet 
and quark jet) and through the boson-gluon fusion graph (where a gluon from the 
hadron splits into a qq pair, one of which scatters off a photon from the incoming 
lepton, to produce a quark and antiquark jet). The ~ 100 < Q 2 <~ 10 4 GeV, 
rapidity, mass, and x distributions of these dijet cross sections are generally well 
described by NLO QCD 0. 

The NLO dijet predictions, which include up to three emitted partons, represent 
LO predictions for three jet production. These LO, three jet predictions are in 
reasonable agreement with HERA three jet cross sections as a function of Q 2 and 
M3j et E3. There does seem to be a discrepancy between data and theory for the ratio 
of three to two jet cross sections. At Q 2 below 1000 GeV 2 the experimental ratio 
is below the LO prediction 113 . However, parton shower simulations show excellent 
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agreement with the data, this is typically a signal that higher order perturbative 
QCD predictions will show improved agreement. The ratio of three jet to two jet 



cross sections has also been measure 
the Et of all jets in the event EEt ^ 



at the Tevatron as a function of the sum of 
The ratio rises rapidly with 'EEt but levels 
off at 0.7 above 200 GeV. The data is fairly well described by LO QCD, but as 
expected is sensitive to the choice of renormalization scale. 

Recently a theoretical milestone was reached with the completion of NLO three 
jet predictions c3. At NLO three jet predictions includes three and four parton 
final states. Figure [| shows a Dalitz lego plot calculated with NLO QCD , where 
Xi = 2Ei/M%j and the indices i = 3,4,5 represent the final state jets ordered in 
energy, Ei is the energy of the three jets, and M$j the invariant mass of the jets. 
The pronounced peaking a x 4 = xz ~ 1 indicates that the third jet is of low energy 
or that the events are primarily two jet in nature. The CDF collaboration has made 
preliminary measurements of the Dalitz distributions and reports good agreement, 
however uncertainty analysis are pending Ea. 

Multijet production also offers an opportunity to search for new dynamics. Al- 
^ 10 ~ 
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Figure 3: Neutral current cross sections for e + p and e p scattering compared to 
Standard Model expectations. 




Figure 4: A Dalitz plot of a NLO three jet cross section. 
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most all QCD data are well described by DGLAP evolution which can be char- 
acterized by partonic emission strongly ordered in Q 2 . The familiar expansion 
a ~ (a s lnQ 2 ) n describes a remarkable range of x and Q 2 . (See the discussion 
of Structure Functions, for example.) However, at fixed Q 2 arid small x, Balitsky- 
Fadin-Kuraev-Lipatov or BFKL evolution may be appropriate i3. Here 1/x becomes 
large and cross sections are of the form a ~ (a s ln(l/x)) n where x = Q 2 /s and s 
is the center-of-mass energy. Terms in a s ln(l/x) correspond to gluon emission 
strongly ordered in x. 

At the Tevatron, for events characterized by two jets widely_separated in A?;, 
the term a 8 ln(l/x) is proportional to a s Ar/. BFKL resummationE3 of these leading 
log terms yields an exponential increase of a with Ar/, a ~ e A '', this approxi- 
mation is valid only for A?y > 2. Crudely, a gluon ladder between the outgoing 
partons increases the cross section. However, the exponential growth is moderated 
by momentum conservation (equivalently the pdf's). Sensitivity to the exponen- 
tial increase can be found in the ratio of cross sections at different center-of-mass 
energies but identical X\ and xi soJiiat the pdf damping cancels. This ratio has 
been measured at Tevatron energies 2 -!! of 1800 and 630 GeV. The data for Ar\ > 2 
is well above and inconsistent with exact LO and HERWIG predictions. (Herwig 
is a LO simulation which includes fragmentation and hadronization) . The BFKL 
prediction best describes the data but is also 3er below the observation. 

HERA is a natural site for BFKL searches since x values near 10 -5 are accessible. 
For example, inclusive tt° production for 10~ 4 < x < 10~ 2 and 2 < Q 2 < 70 GeV 2 
is well described by LO BFKL predictions but poorly described by LO QCD E3. 
Likewise forward jet production or production in thcJheam direction, 1.5 < r\ < 2.8 
and 10 -4 < x < 10~ 2 is not described by NLO QCD t3. There is also a pronounced 
dependency on the choice of renormalization scale. This suggests that the calcula- 
tions are incomplete and something more is needed, perhaps BFKL improvements, 
a better description of the photon structure function, or next-to-next-to-leading 
order terms. 

BFKL dynamics are also accessible at LEP through 7*7* scattering where each 
incoming lepton radiates a virtual photon. At LO the photons scattering off a 
qq pair. At higher orders each photon emits a qq pair and the two pairs interact 
through gluon exchange. The L3 experiment at LEP has measured a rise in (7(7*7*) 
as a function of the scattering center-of-mass c3. This increased cross section can be 
attributed to the the phase space opened up by emission of additional gluons from 
the exchanged gluon. The rise is not described by LO QCD but is described by a 
BFKL prediction. There are abundant hints for the existence of BFKL dynamics; 
however, all searches need improved higher-order calculations, both DGLAP and 
BFKL, before the origins of these hints are clear. Moreover, DGLAP has been 
impressively successful over almost all x and Q 2 . 

4. The Structure of the Proton 



The proton pdf's are derived from global fits to two general types of data: structure 
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functions from e, fi and v scattering from nuclear targets and from exclusive state 
production such as Drell-Yan, charm, W, jet, and photon production. Without 
going into details, trial pdf's are evolved form a starting Q a , convoluted with QCD 
calculations of the hard scatter, and fit to the observed data. These pdf's are 
assumed to be universal in that a single unique pdf is appropriate for all reactions. 
Each reaction does, however, provide access to a particular type of pdf. For instance, 
W production in pp scattering directly accesses the u and d quark pdfs. Generally, 
quark functions are heavily constrained by DIS scattering for x < 0.9. Until recently 
there has been very little to constrain gluons above x > 0.2 (This freedom 
permitted the formulation of pdf's similar to CTEQ4HJ used in Fig. |[) Recent 
developments in pdf derivations include an increased accuracy due to new types of 
input data, replacement of prompt photon data with jet data, and concerted efforts 
to improve the treatment of uncertainties. For a summary consult the references Ej. 

The differential cross section d 2 a NC /dxdQ 2 for NC scattering, where x is the 
momentum fraction of the struck quark, can be expressed at leading order in terms of 
the structure functions F2(x, Q 2 ) and i<3(x, Q 2 )- F 2 is proportional to the sum of the 
valence and sea quark pdf's, q(x, Q 2 ) + q(x, Q 2 ). F 3 is proportional to the difference 
of the valence and sea quark pdf's, q — q. Figure I shows F 2 for 0.000032 < x < 0.65 
and 1 < Q 2 <~ 10 4 GeV 2 as measured by the HI, SLAC, NMC and BCDMS 
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Figure 5: A compilation of structure functions over a wide range of x and Q 2 . 
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collaborations. (To aid the eye the data at each value of x is offset by an arbitrary 
constant.) The HI data is recent, and ZEUS has made similar measurements El. 
The curves correspond to a NLO QCD fit and include both e + p and e~p data . 
The astonishing agreement includes both the large Q 2 perturbative region and also 
the small, presumably, nonperturbative Q 2 regions. The dependence of F 2 on Q 2 
at low-x region is an example of QCD scaling violations. Although not shown, the 
availability of both e + p and e~p NC data also permits measurement of -F3, these 
results are in agreement with NLO QCD El 

Charge current reactions couple directly to the quark distributions. Specifi- 
cally, d 2 a cc \e~ p) / dxdQ 2 is proportional to the u and c quark distributions, and 
d 2 a cc (e + p)/dxdQ 2 is proportional to the d quark. Figure ^| shows a reduced CC 
cross section for various Q 2 bins as a function of x for e~p data El The uppermost 
curve in each panel represents the total cross section as predicted by NLO QCD us- 
ing the CTEQ5D pdf. Agreement is good at all x and Q. The lower curves indicate 
the contribution from each type of quark. As expected the u quark contribution is 
dominant. 

Additional constraints on the quark distributions below x ~ 0.7 can be de- 
rived from fixed target and pp scattering. A few recent examples are listed here. 
The CCFR collaboration has re-analyzed its v and v DIS data and finds F2 for 
0.015 < x < 0.175 and 1 < Q 2 < 100 GeV 2 well described by updated NLO QCD 
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Figure 6: The CC cross-section at various x and Q 2 with predictions for the various 
quark contributions. 




A QCD Survey: < Q 2 < 10 5 GeV 2 9 

predictions El. CCFR has also made the first measurement of F£ — F% . The rapidity 
distributions of the final state leptons in the Drell-Yan process pp — » 7* /Z — » e + e~ 
is sensitive to the valence quark distributions. For example, if the two initial par- 
tons are of unequal momentum the virtual 7*/Z will be boosted in the direction of 
the more energetic parton and the decay electrons will likewise be boosted. Anal- 
ysis of the rapidity distributions as recently measured by_the CDF collaboration 
constrains the u and d distributions for 0.05 < x < 0.6 l3. Recent results also 
pertain to the sea-quark distributions, the cross sections for [i + \i~ production from 
pp and pd scattering as measured by the E866 collaboration constrains the d — u 
sea distribution between x ~ 0.03 and ~ 0.3 Ej. 

As mentioned earlier, the large-x gluon distributions are not well known. Until 
recently, a favored process for constraining g(x) in the global fits was direct photon 
production. The associated Compton graph samples g(x) directly. However, this 
process has fallen out of favor since NLO calculations cannot describe the world 
data 123. For nearly every photon cross section measurement, the data exceeds the- 
ory at the low end of that specific measurement's x range. This discrepancy has been 
attributed to initial-state soft radiation which manifests itself as transverse momen- 
tum or kr- If not interpreted correctly, this kx or "kick" adds to the measured x of 
the photon which creates the observed excess of data over theory. Evidence for kx 
comes from events where object pairs such as pions, diphotons, dimuons, and dijets 
are produced. The vector sum of the pair momentum is a direct measure of the 
parton scatter transverse kick. Indirect evidence comes from the observation that 
augmentation of NLO predictions with a phenomcnological addition of fc-r dramat- 
ically improves agreement between data and theory. Also, the importance of soft 
gluon resummation for Drell-Yan, diphoton, W, and Z production was recognized 
some time ago. Similar calculations for direct photon production are underway and 
already show much better agreement between data and theory. Perhaps with time 
and higher orders photon data will return to the pdf fits. 

In the meantime, the central inclusive jet cross sections constrain g{x) for ~ 
0.05 < x <^ 0.5. Inclusion of the full inclusive cross section for pp scattering over 
7/| < 3.0 will cover x <~ 0.8 □. Even more can be gained from triple differential 
cross sections, d 3 a/dETdrjidr]2 where Et corresponds to the leading jet in an event 
and the subscripts identify the leading two jets u. For a dijet event if both jets 
are at r\ = with Et = 180 GeV then x\ = x-i — 0.2. However, for r\ = and 
Et = 90 GeV with the second jet at 77 = 2, x\ = 0.2 and x% =~ 0.9. Since the event 
structure is boosted forward a much greater fraction of the initial hadron momentum 
is required. Both Tevatroii^experiments are encouraged to publish their preliminary 
differential cross sections EltJ since NLO predictions exists and the data could be 
incorporated into the pdf's. Dijet production at HERA, through the boson fusion 
graph, also directly measures g(x) < 0.1 Ej. Although the data doesn't greatly alter 
the pdfs's relative to the F2 data it does add stability to the fits and reduce errors. 

Until recently the uncertainties of observables due to pdf's were estimated merely 
by using a menu of pdf's. However, pdf's should more properly handle statistical 
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and systematic data uncertainties, theoretical uncertainties, and parameterizations. 
In the past yea r or two, there has been very great progress dealing with these 
uncertainties EjOE3. For instance, the plot below by Botje shows the allowed 
variations in the gluon and quark distributions due to statistical uncertainties, input 
choices, analysis cuts, and renormalization choices assuming all uncertainties are 
symmetric and Gaussian in nature O. Notice the unconstrained nature of the gluon 
distributions, Tevatron jet results were not included in this calculation. An even 
more general approach by Giclc and Keller incorporates probability distributions to 
derive the pdf 's E3. These authors and others are providing tools for incorporating 
all sources of pdf uncertainties. 



5. Status of a s Measurements 

An enormous body of research has been dedicated to the study of the strong coupling 
constant, a s since it is the only free parameter of QCD and must be determined 
experimentally. Of equal interest is the dependence of the coupling constant on 
momentum transfer Q 2 , a s (Q 2 ) = 12tt/ (33 — rif)log(Q 2 /A 2 ) where rif is the number 
of quark flavors and A is experimentally determined. Notice a s decreases or "runs" 
with momentum transfer. This is, in fact, the basis for the perturbative NLO QCD 
calculations described earlier. 

The strong coupling constant can be derived in a myriad of ways and at all Q 2 , 
from absolute decay rates of the Z boson and r lepton, energy levels of bound heavy 
quarks, jet event shapes, jet production rates and angular distributions, and scaling 
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Figure 7: Allowed variations of the gluon and quark pdf's assuming symmetric 
Gaussian uncertainties associated with the global fit. 
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violations in deep inelastic scattering. Notable new measurements of a s originate 
from all three high energy colliders. In particular, the LEP collaborations recently 
ran at center-of-mass energies from 192 to 202 GeV and from event shape and 
jet rate measurements determined an average a s (198GeV") = 0.109 ± 0.001 ± 0.005 
where the first error is statistical and the second theoretical til. A uniform re- 
analysis of e + e~ event shapes, jet rates, and multiplicities from the JADE and 
LEP collaborations from center-of-mass energies between 20 and 200 GeV show the 
strong coupling constant obviously running with a value of a s (Mz) = 0.1208 ± 
0.0006 ± 0.0048 E3. Finally, comparison of NLO inclusive jet and dijet cross section 
predictions to experimental results from HERA and the Tevatron consistently yield 
a s (Mz) measurements near 0.12 US The Tevatron inclusive jet measurement 
beautifully demonstrates the running of a s over 10 3 < Q 2 < 2 x 10 5 GeV 2 . 

An excellent review by Bethke includes a comprehensive compilation of the many 
derivations of a s complete with a description and evaluation of each measurement 
l3. The current world average for the coupling constant, a s (Mz) = 0.1184±0.0031, 
is based on six measurements for which next-to-next-to-leading order predictions 
exist. The average is known to better than 3% and is unchanged since 1997 0. 

6. Conclusions 

QCD, and in particular perturbative QCD, describes the strong interaction over an 
impressive kinematic range. At the highest Q 2 accessible, recent data and analyses 
reveal few, if any, excursions from the Standard Model. At intermediate Q 2 , jet and 
multi-jet studies have permitted a closer study of higher order process and searches 
for BFKL dynamics. These searches are ambiguous. The proton structure functions 
are under intense study. Although the quark components are well constrained, the 
gluon component is not well measured for x > 0.2. Recent jet and, if rehabilitated, 
photon production data should improve this situation. There has also been great 
progress on pdf error analysis. Data from all values of Q 2 have contributed to a 
precise 3% measurement of the strong coupling constant. 
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